Quantum localization in bilayer Heisenberg antiferromagnets with site dilution 
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The field-induced antiferromagnetic ordering in systems of weakly coupled S = 1/2 dimers at 
zero temperature can be described as a Bose-Einstein condensation of triplet quasiparticles (singlet 
quasiholes) in the ground state. For the case of a Heisenberg bilayer, it is here shown how the above 
picture is altered in the presence of site dilution of the magnetic lattice. Geometric randomness 
leads to quantum localization of the quasiparticles/quasiholes and to an extended Bose-glass phase 
in a realistic disordered model. This localization phenomenon drives the system towards a quantum- 
disordered phase well before the classical geometric percolation threshold is reached. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Cx, 64.60.Ak 



The introduction of randomess in strongly correlated 
fermionic and bosonic systems on a lattice offers a large 
showcase of novel quantum behavior In quantum 

spin systems vifith an ordered ground state, the generic 
effect of geometric disorder is an enhancement of local 
quantum fluctuations through the reduction of the con- 
nectivity of the lattice 0|, which may lead to quantum 
disordered phases 0|. On the other hand, in spin sys- 
tems with a gapped quantum-disordered ground state, 
geometric randomness can induce local magnetic mo- 
ments giving rise to long-range order through an order- 
by-disorder mechanism In the case of superfluid in- 
teracting bosons, the presence of an increasing amount of 
disorder can destroy superfluidity through quantum lo- 
calization, inducing a disordered, yet gapless Bose-glass 
phase with a finite compressibility . 

Recently, the strong connection between lattice bosons 
and quantum spin models has been pointed out in the 
context of unfrustrated 5=1/2 weakly-coupled dimer 
systems. For sufficiently weak inter-dimer couplings, the 
spin system has a singlet ground state with a finite gap 
to triplet excitations. The application of a uniform mag- 
netic field leads to a lowering of the energy of the triplet 
excitations aligned with the field, until, at a lower crit- 
ical field hf^, the ground state develops a finite mag- 
netization. Within a bosonic mapping that neglects 
the other two triplets, bosonic hardcore triplet quasi- 
particles (TQPs) populate the ground state forming a 
bosonic condensate |G|. Off-diagonal long-range order 
corresponds to the appearence of a finite staggered mag- 
netization in the plane perpendicular to the field. When 
the magnetization is half-saturated, the relevant degrees 
of freedom of the condensate become hardcore singlet 
quasi-holes (SQHs) in the triplet sea, which disappear 

completely at an upper critical field /ij;2^ fully polarizing 
the system. For h < h^^j^ {h > h^!^) the system is there- 
fore a band insulator, with an empty (full) band of TQPs. 
The validity of the description of the ordered phase as a 
Bose-Einstein condensate has been extensively investi- 
gated both theoretically and experimentally 0; S ■ 

By virtue of the bosonic picture, it is intriguing to 



envision the possibility of a Bose-glass phase for weakly 
coupled dimer systems in a magnetic field upon introduc- 
tion of lattice disorder This opens up the appealing 
perspective of unambiguously realizing such a phase in a 
quantum spin system, profiting of the high level of con- 
trol with which disorder can be introduced in a magnetic 
lattice [iJl . The focus of this Letter is the extent and the 
major signatures of such a phase in a realistic model. 

In what follows we investigate the S — 1/2 Heisenberg 
antifcrromagnet in a magnetic field and with site dilution, 
whose Hamiltonian reads: 

(ij) a=l,2 

+ ■ Si^2 ~ gfJ^sH'^ei^aSl^ . (1) 

Here the index i runs over the sites of a square lattice, 
(ij) are pairs of nearest neighbors on the square lattice, 
and a is the layer index. J is the mterlayer coupling and 
J' the mtralayer one. The variables ei^a take the values 
or 1 with probability p and 1—p respectively, p being the 
concentration of non- magnetic sites. Hereafter we will ex- 
press the field in reduced units h — gfiBH/ J . We have in- 
vestigated the above Hamiltonian making use of Stochas- 
tic Series Expansion quantum Monte Carlo based on the 
directed-loop algorithm and considering L x L x 2 
lattices with L up to 40. The T — limit is reached 
through a /3-doubling approach Q. Disorder averaging 
is performed by using typically 200 different disorder re- 
alizations. 

In the clean limit {p — 0) and at zero field, the Hamil- 
tonian of Eq.ffl) is in a gapped dimer-singlet phase for 
J /J' ^ 2.5 [13. Hereafter we will specialize to the 
case of a bilayer with J/ J' = 4, namely well inside 
the dimer-singlet phase. As discussed in the introduc- 
tion, the clean system with large J/ J' ratio can be 
well described by an effective model of hardcore TQPs 
which hop from dimer to dimer with amplitude J'/2, 
nearest-neighbor repulsion J'/2, and chemical potential 
J{h - 1) controlled by the field 0. Increasing the 
field leads to the closing the triplet gap, accompanied 



2 



by the condensation of TQPs. We have investigated this 
scenario numerically, by studying the field evolution of 
the uniform magnetization mf^ — (Sf), the uniform sus- 
ceptibility Xu = dm^/dh, the staggered magnetization 



j), and the 



superfluid density T = T/2 {W^ + W^), where 

is the worldline winding number in the x{y) direction. 

As shown in Fig. ^ a field-induced ordered phase with 



rrig, T 7^ is clearly observed between hf^ = 0.47(1) and 
^c2 = 2, in between a dimer-singlet phase {h < ft.^'j^) and 
a fully polarized phase {h < hf^). 
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FIG. 1: Field dependence of the uniform magnetization m^, 
transverse staggered magnetization mj, and superfluid den- 
sity T at r = for the bilayer antiferromagnet with J/ J' = 4 
in the clean limit (p = 0). 

We now turn to the case of the doped system, p ^ 0. 
From the purely geometric point of view, the lattice un- 
dergoes a percolation transition at p* = 0.5244(2) which 
we estimated through classical Monte Carlo calculations 
based on a recently proposed algorithm . 

The magnetic behavior in turn displays a rich se- 
quence of phases well below the percolation transition. 
Fig. [21shows the zero-temperature phase diagram, where 
boundary lines between ordered and disordered phases 
have been estimated through the linear scaling of the 
disorder-averaged correlation length, (j^^'^yy) ~ L. 

In what follows, we will discuss the phase diagram 
making use of both the bosonic and the magnetic lan- 
guage. Strictly speaking, the conventional bosonic map- 
ping onto hardcore TQPs breaks down in presence of im- 
purities, since some of the spins become unpaired^ as they 
miss their partner on the neighboring layer. A bosonic 
mapping is however still possible in which bosons corre- 
spond to I tt) states on intact dimers and to | t) states 
on unpaired spins. Site dilution of the spin model re- 
flects in correlated diagonal and off-diagonal disorder of 
the bosonic model [l6|. 

In the dimer-singlet phase the introduction of vacan- 
cies in the magnetic lattice leads to the appearence of 
local 5 = 1/2 moments. It is straightforward to iden- 
tify these moments with unpaired spins, but in fact they 




FIG. 2: Ground-state phase diagram of the site-diluted bi- 
layer Heisenberg antiferromagnet with J / J' = 4. 




FIG. 3: (a)-(b): Real-space images of the dimer magnetiza- 
tion rui — {Sli + Sf^2) on intact dimers in a 40x40x2 bilayer 
with J /J' = 4, dilution p = 0.1 and at inverse temperature 
PJ = 256, for h = 0.56 (a) and h = 0.6 (b). The radius of the 
dots is proportional to the dimer magnetization. The mag- 
netization of unpaired spins is omitted for clarity. The most 
visible localized states are highlighted in (a), while the back- 
bone of the percolating magnetized network is highlighted in 
(b). (c): Superfluid density as a function of the held for the 
specific sample considered. 

occupy a larger volume ^ where £o is the correla- 
tion length in the clean limit p = |l7l |. thus spread- 
ing also over intact dimers. Those moments are coupled 
through an effective unfrustrated antiferromagnetic net- 
work, which can sustain long-range order at T = 0. This 
gives rise to an order-by-disorder phenomenon, which 
is expected to persist for any doping concentration up 
to percolation p*. The network of local moments has 
a broad distribution of effective couplings, which scale 
exponentially with the inter- moment distance r Il8l |. 
JcS ~ exp[— r/^o]- The application of a moderate field 
(h < h^^^) can easily destroy the antiferromagnetic order 
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by polarizing a majority of the effective magnetic mo- 
ments |l9j , bringing the system into a partially polarized 
quantum- disordered phase. In particular, those dimers 
which are immediately close to an unpaired spin (roughly 
within a distance of ^o) are also partially polarized, even 
if the field is lower than the dimer gap and the bilayer 
gap. In the bosonic picture, this means that TQPs are 
present in the ground state, localized around unpaired 
spins and only weakly penetrating clean regions of intact 
dimers, which mostly remain unmagnetized. They form 
a "gossamer" superfluid network with exponentially sup- 
pressed links for small fields, but for larger particle num- 
bers this network turns into an insulating one due to the 
complete filling of some of its islands. 

Upon increasing the field beyond h^?^ , the picture 
of TQPs exponentially localized around vacancies start 
to change substantially. Continuous regions of intact 
dimers, which occur with probability exponentially small 
in the size of the region, begin to respond to a field simi- 
larly to what would happen in the clean case. Some of the 
intact dimers with a lower local coordination, nonethe- 
less, have a higher local gap which a field h w h^^j^ is not 
able to close, so that they remain non-magnetized up to 
higher fields. In the bosonic picture, this means that 
TQPs begin to appear in the bulk of the clean regions, 
but they are exponentially localized due to disorder. This 
picture clearly corresponds to that of a Bose-glass phase 
Such a phase extends up to a critical field hd at 
which the increased population of bosons, together with 
the hardcore repulsive interaction, leads to a localization- 
delocalization transition and to superfluidity, as shown in 
Fig. 131 There the density distribution of the TQPs is im- 
aged by the dimer magnetization rrii — {S^i + 3^2)- 
the figure shows, superfluidity occurs via quantum per- 
colation of the collective bosonic state throughout the 
lattice, while single particle states would all be exponen- 
tially localized in 2D. 

From the numerical data, we observe that the Bose- 
glass phase extends over a field region that increases with 
the dilution p. In particular, beyond an upper critical di- 
lution, pc ~ 0.36, no field-induced order survives. This 
happens well below the lattice percolation threshold p* . 
Therefore there is quite a large region of doping com- 
pletely dominated by quantum disorder. 

As the field is increased even further for p < pc, sin- 
glet quasi-holes (SQHs) become the relevant degrees of 
freedom. They undergo a similar superfluid-to-Bose-glass 
transition at a field hc2 without major alterations with 
respect to the TQP case. In particular the magnetization 
does not reach full saturation until the clean upper crit- 
ical field h^^2 y given that (exponentially rare) large clean 
areas cannot be fully polarized before h^^2 is reached. 

A typical field scan through the phase diagram is 
shown in Fig. ^for a dilution p — 0.2. The field is seen 
to quench the order-by-disorder spontaneous transverse 
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FIG. 4: Zero-temperature field scan in the site-diluted bilayer 
Heisenberg antiferromagnet with J/ J' — 4 and p = 0.2. Inset: 
scaling of the order parameter in the different phases of 
the system; from top to bottom h = 0.75 (XY ordered), 0.05 
(order- by-disorder), 1.85 (high-field Bose glass), 0.35 (quan- 
tum disordered), 0.6 (low-field Bose glass). 

magnetization and to drive the system towards a partially 
polarized state without fully polarizing the local free mo- 
ments, whose magnetization continues to grow even be- 
yond the transition. Eventually for large enough fields 
all the local moments get fully polarized by the field, 
leading to a magnetization plateau at to = pS = p/2 
[20| . In the Bose-glass phase the uniform magnetiza- 
tion is seen to change concavity and to start to grow 
slowly, as the clean regions are gradually magnetized. 
A simple analysis can be made based on a "local-gap" 
model, for which the overall magnetization is the sum 
of the magnetizations of uncorrelated clusters of intact 
dimers with different local gaps jl^. The result is that 
the growth over the plateau value is exponentially ac- 



tivated. 



pS 



exp 



where the 



exponential behavior comes from the exponential tail of 
the size distribution of clean regions, and it is a direct 
evidence of the dominant behavior of rare events in the 
Bose-glass phase. Similarly, an exponential saturation is 



observed for h 



,(0) 



t^2 in the upper Bose-glass region. 

Let us now turn to the finite-temperature signa- 
tures of each phase. The zero-temperature XY ordered 
phase turns into a quasi-long-range ordered phase at fi- 
nite T up to a Berezinskii-Kosterlitz-Thouless (BKT) 
transition temperature T^^^, at which the superfluid 
density vanishes and the transverse structure factor 

S{n,n) = (l/A^)E.,,./5(-l)'+-''+"+^(C^^S'^) be- 
comes finite. We estimate the transition temperature 
Tj3j^^ through BKT scaling of the superfluid density and 
of the transverse structure factor [2l|. The transition 
temperature as a function of the fleld is shown in Fig. 
Due to disorder and quantum fluctuations, T^j^^ is 
almost two order of magnitudes lower than the interlayer 
exchange coupling J, and its slightly asymmetric field 
dependence mimics that of the T = order parameter. 
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FIG. 5: Finite-temperature phase diagram for the site-diluted 
bilayer Heisenberg antiferromagnet with J /J' =4 and dilu- 
tion p — 0.2. In the inset: temperature dependence of the 
transverse structure factor for the same model with L = 32. 



Turning on the temperature from a T — disor- 
dered phase, on the other hand, we observe interesting 
differences between the quantum disordered phase for 
h < h'^i and the Bose-glass phase for h'li < h < hd 
and hc2 < h < /i°2- In the clean Umit p = 0, the gapped 
disordered phases for h < h^^ and h > /i°2 show the com- 
mon feature of thermal activation, with transverse struc- 
ture factor S'(7r,7r), transverse correlation length ^^^(yy\ 
and uniform magnetization increasing as T grows 
above 0. Such a behavior is easily understood: upon 
increasing the temperature, the system is able to explore 
more correlated states with finite magnetization sitting 
beyond the energy gap. In the presence of site dilution, 
the temperature activation of correlations turns into a 
slow decrease as T increases in the low-field regime dom- 
inated by the unpaired spins (see inset of Fig. |SJ). For 
high enough fields, and in particular in the Bose-glass 
phase, the temperature activation of correlations is re- 
stored up to a temperature Tg for ^(Tr, tt) and for ^ 
at which both quantities hit a maximum, showing then a 
crossover to a decreasing dependence typical of a conven- 
tional paramagnetic phase. In contrast to the clean disor- 
dered phases, the Bose-glass phase is gapless. Nonethe- 
less some clusters of intact dimers resist being magne- 
tized for h^^^ < h < hci or they are fully polarized al- 
ready for hc2 < h < . This means that the observed 
temperature activation in that field range can be asso- 
ciated with those portions of the system visiting states 
with higher transverse correlations, linked to the pres- 
ence of bosonic TQPs/SQHs. In the bosonic language, 
this can be suggestively pictured as a partial delocal- 
ization of the trapped quasiparticles into those regions, 
assisted by thermal fluctuations. 

Our findings are directly relevant to experimental in- 
vestigations of site-diluted weakly coupled dimer systems 
in a magnetic field. The particular case of an Heisen- 
berg antiferromagnet on coupled bilaycrs is realized by 



BaCuSi2 06 [g, in which doping of the magnetic Cu^"*" 
ions with non magnetic Zn^+ and Mg^+ can lead to site 
dilution of the magnetic lattice. The picture of quantum 
localization of bosonic quasiparticles applies nonetheless 
to other spin gap systems with different geometries, such 
as Sr2Cu(B03)2 ^ or Tl(K)CuCl3 0. The quanti- 
ties we have investigated clearly show a characteristic se- 
quence of phases, with disorder suppression of magnetism 
and onset of a novel Bose-glass phase, and they are di- 
rectly accessible to magnetometry and neutron scattering 
experiments. 
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